Three-Particle Correlations in Simple Liquids 
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We use video microscopy to follow the phase-space trajectory of a two-dimensional colloidal model 
liquid and calculate three-point correlation functions from the measured particle configurations. 
Approaching the fluid-solid transition by increasing the strength of the pair-interaction potential, 
one observes the gradual formation of a crystal-like local order due to triplet correlations, while being 
still deep inside the fluid phase. Furthermore, we show that in a strongly interacting system the Born- 
Green equation can be satisfied only with the full triplet correlation function but not with three- 
body distribution functions obtained from superposing pair-correlations (Kirkwood superposition 
approximation) . 
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Our current understanding of the structure of simple 
fluids is based on the n-body distribution functions g^ n \ 
measuring the probability density of finding two, three, 
and more particles at specified positions in space. When 
the total potential energy of a liquid is given by a sum of 
pair-potentials, all of its thermodynamic properties can 
be calculated by means of the pair-correlation function 
g(r) = g( 2 '(r) and its density (p) and temperature (T) 
derivatives. However, the latter two quantities, dg(r)/dp 
and dg(r)/dT, explicitly depend on the triplet correla- 
tion function, even if the particle interactions are only 
pair- wise additive Q] . Explicit knowledge of triplet corre- 
lations is also required in perturbation theories for static 
fluid properties |2J, in theories of transport properties 
3, of solvent reorganization processes around solutes 
ij, of systems under shear-flow ||, but also to under- 
stand the structural properties of a 2D amorphous sys- 
tem 0. Most of our knowledge on triplet correlations 
come from computer simulation studies of hard-sphere- 
fluids RJ, Lennard- Jones fluids 0, E| and electrolyte sys- 
tems jlflj. In the overwhelming majority, these papers 
are concerned with testing Kirkwood's superposition ap- 
proximation (KSA) for the triplet distribution func- 
tion. By contrast, semi-analytical theories for g^ be- 
yond the KSA are rather rare 0,0. However, despite 
the long-standing theoretical interest in its properties, it 
has never been possible to measure three-particle corre- 
lations directly. Indirect ways to identify higher-order 
correlations in scattering data have been suggested for 
instance in [l3| . An alternative, but also indirect way 
to obtain experimental information on g^ is based on 
the relationship between the isothermal pressure deriva- 
tive of the fluid structure factor dS(q) / dP and the triplet 
distribution function 0, a relationship which has been 
systematically exploited by Egelstaff and co-workers in 
rare-gas systems |l4J. The present Letter reports on the 
first direct measurement of g^ in a two-dimensional col- 
loidal model liquid with well-defined pair-interaction po- 
tentials. 

The preparation of the samples and the experiments 
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FIG. 1: A typical image (500 x 380 fim) of our two- 
dimensional colloidal model system with paramagnetic col- 
loids of d = 4.7/im diameter. The particles interact via a 
potential ~ Y/r 3 in which the interaction strength F can be 
conveniently varied through the external magnetic field. 



have been carried out as described in |T^: Spherical 
colloids (diameter d = A.l /im) are confined by gravity 
to a water/air interface whose flatness can be controlled 
within less than a micron. The field of view has a size 
of 520 x 440 [vm containing typically about 10 3 particles. 
The particles are super-paramagnetic due to Fe2 03 dop- 
ing. A magnetic field B applied perpendicular to the 
air/water interface induces in each particle a magnetic 
moment M — \B which leads to a repulsive dipole-dipole 
pair-interaction energy of f3u(r) = r/(y / 7rpr) 3 with the 
interaction strength given by T — /3(/! /47r)(x-B) 2 (7r/9) 3 / 2 
(/3 = 1/kT the inverse temperature, x the susceptibility) . 
This is the only relevant contribution to the interparticle- 
potential which is hence conveniently and reversibly ad- 
justable by the external field B. A typical snapshot of our 
system is given in Fig. V is the only parameter deter- 
mining the phase-behavior of the system: for Y < 57 the 
system is liquid, for Y > 60 it is solid, and in between, 
i.e. for 57 < Y < 60, it shows a hexatic phase [l5j |. 



We here analyze three different T(r = 4, 14,46), where 
the system is deep in the liquid phase, and use for each 
(well equilibrated) system about 200 statistically inde- 
pendent configurations with approximately 500 particles, 
recorded using digital video-microscopy with subsequent 
image-processing on the computer. From the measured 
particle configurations, is obtained by computing the 
average count per configuration of a particular kind of 
triplet, divided then by the appropriate normalizing fac- 
tor. Details of this calculation will be given elsewhere 

M 

Triplet correlations can be characterized by the ra- 
tio between the full triplet distribution function g^ 

(3) 

and its approximated form based on the KSA g SA = 
g(ri)g(r2)g(r^). This ratio is given by what is 
called the triplet correlation function, denoted here by 
G(ri, r2, r3). Thus, g^ — gf^G. All pair correlations in 
g^ are included in gg A , while the extent of the intrinsic 
correlations due to the simultaneous presence of a triplet 
of particles at positions ri, T"2, r% is quantified through 
the function G, which thus defines the local structure of 
the fluid beyond that expressed by the pair-correlation 
functions. Introducing /3w^ m ^ = — \o.g^ m \ g^ = gf A G 
transforms into 

w (3) (run, r 3 ) = w (2) (n) + ™ (2) (ra) + w (2) (r 3 ) - In G/(3 

(1) 

with and being the two- and three-particle po- 
tential of mean force, respectively. The equation shows 
that — In G = flAw^ plays the role of a three-body po- 
tential, measuring the (extra correlation) energy of three 
correlated particles relative to the energy of superposed 
correlated pairs of particles. 

In an homogeneous, isotropic system, g^ depends 
on only three independent variables, chosen here to be 
r = |ri - r 2 |, s = \r 2 - r 3 | and t = \r x - r 3 \. Fig. @ 
shows three-particle distribution functions in the equi- 
lateral triangle geometry for all three T's considered. To 
allow comparison with the radial distribution function 
g(r) we have taken the cubic root y/ g^(r, r, r) so that 

5si( r ' r ' r ) = ff( r )- I* i s evident that the KSA, while 
working satisfactorily at low T, fails to reproduce the fine- 
structure of the triplet distribution function at higher 
values of T. Obviously, correlations beyond the level of 
pair-correlations become important at higher T. To visu- 
alize g^ in two dimensions, r is fixed in the following to 
(2) 

the distance r max where g[r) has its first peak (« 1/^/p), 
so that g^H r< max, s,t) varies just with s = s(x,y) and 
t = t(x,y) and can thus be plotted in the (2;, y)-plane 
in form of a contour-plot. This is done in Fig. for 
the r = 46 measurement. We show in the left half of 

(3") 

the figure g SA and contrast it to the full three-particle 
distribution function g^ 3 \ plotted in the right half of 
the figure, approaches g(rmax) for large values of 

x 2 + y 2 . To keep the figure as clear as possible, we plot- 
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FIG. 2: Triplet distribution functions </ 3 ' as a function of 

the side length of an equilateral triangle, as computed from 

(3) 

measured particle configurations for different T. g SA is the 
triplet distribution function in the Kirkwood superposition 

approximation which on taking the cubic root, y g t g\( r ' r > r )> 
becomes the radial distribution function g(r). 



ted just those parts of g^ - g(r ( m ' ax ) and g ( SA - g(r ( m ' ax ) 
that are larger than zero. The stripes that can be seen 
especially close to the x-axis result from the transforma- 
tion g^(r max ,s,t) to g( 3 Hr max ,x,y) and appear due to 
limited statistics. A hexagonal lattice with a lattice con- 
stant a = rma X is superposed. Both distributions g^ and 

(3) 

9sa revea l that the neighbors of the two central particles 
have positions which show a certain correspondence to 
the crystalline lattice points. However, while the banana- 

(3) 

like structure of g s { reflects just the coordination shells 
of the lattice, the full distribution function g^ shows 
a well-developed, angular dependent substructure, with 
individual peaks for every lattice point in the first coor- 
dination shell. Fig. (@J shows g^ and g^ A of Fig. © 

(2) (21 

along the line (r = r^Lx, s = r; n ' ax ,t = t((f>)), which is a 
circle of radius r^ax around the right particle in Fig. 
passing through all lattice points of the particle's first 
coordination shell (arrows in Fig. Q mark positions of 
lattice points). It can be clearly seen that g^ 3 ' develops 
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FIG. 3: Distribution functions g (r = rLL, s(x, y),t(x,y)) 
(right half of the figure) and g^(r = r^l x , s(x,y),t(x,y)) 
(left half of the figure) in the (x, y)- plane (F = 46). Fhe 
missing half of each distribution is just the mirror-image of 
the one actually plotted. The constant g(rmax) is subtracted 
from the distributions and only positive values are plotted 
with a grey-level scheme between white (zero) and black (max. 
value). 




FIG. 4: g (3) and gf^ from Fig. © for fixed values of r = Vmax 
and s — Tmax: as a function of the angle (j> (see inset). Symbols 
(solid lines) for distributions generated from measured (MC- 
simulated) configurations. Also given is the logarithm of the 
triplet correlation function G, which is related to the triplet 
correlation energy Au' 3 ', for V = 4 (dashed line), T = 14 
(dashed-dotted line) and F = 46 (dashed line). 



peaks at the lattice points while g s ^ completely fails to 
reflect the hexagonal structure. Also given is the function 
— In G, i.e. Aw' 3 ' , of eq. Q, now for all three values of T 
studied here. It is evident how Au/ 3 ) gradually forms on 
increasing T, with values up to one kT (in other regions 
of the (r, s, t)-space we find energies as high as 4 kT\). It 
is also seen that the regions of attractive and repulsive 
correlation energies Au/ 3 ) correspond to the correcting 
effect which the function G has on g^\ to ensure that g^ 
adapts locally to the hexagonal symmetry. We conclude 
that it is an effect entirely due to three-particle corre- 
lations, i.e. due to the function G, which is responsible 
for the observed formation of a crystal-like local environ- 
ment around particles well below the freezing transition. 
We also performed Monte-Carlo (MC) simulations using 
the above-given pair-potential 0u(r) = r/(y/irpr) 3 (with 
a better statistic than in the experiment: 500 config- 
urations with 2000 particles, periodic boundary condi- 
tions). The almost perfect agreement between the dis- 
tribution functions based on the MC-data (solid lines in 
Fig. Q) and on the experimental configurations (sym- 
bols in Fig. J3J), demonstrates that our model liquid con- 
sists of particles interacting solely via pair-wise additive 
and precisely known potentials. 

Furthermore, we carried out MC simulations for the 
solid phase (r = 80), starting from a perfect hexago- 
nal lattice, and compared the resulting triplet distribu- 
tion function with the experimental one for the T = 46 



measurement in the liquid phase, see |16(. The distri- 
butions look quite similar: as regards the correlations 
between the central pair and the first coordination shell 
(in a plot like that in Fig. 0), there is hardly any dif- 
ference between the liquid and the solid phase. Pro- 
nounced differences are observable, however, in the sec- 
ond shell: in the liquid phase the next nearest neighbors 
are broadly distributed midway between adjacent lattice 
nodes (see Fig. JjJJ), while the T = 80 distribution corre- 
lates much better with the lattice structure. However, 
even for r = 80 this correspondence is far from per- 
fect; it is well developed along the y-direction, but be- 
comes worse on increasing <f> to 180° where there is still 
an extended smeared-out distribution showing no clear 
preference for certain lattice points. Clearly, approach- 
ing T — > (r — > oo), one will ultimately observe peaks 
in g^ positioned exclusively on the lattice points. We 
should remark that in three dimensions a similar corre- 
spondence between the peaks in g^ and an underlying 
crystal lattice should be much harder to find. In 3D, ev- 
ery triplet of particle lies, of course, also in a plane, and 
can accordingly be plotted as in Fig. J3J . However, then 
there is not one, but a superposition of many possible 
lattice planes that one has to compare this distribution 
with. 

To demonstrate that triplet correlations are significant 
not only locally, but also when integrated over the whole 
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FIG. 5: Fest of Kirkwood's approximation using experimen- 
tally determined three-particle distribution functions (r = 46 
and r = 4). Solid lines for the left hand side of the Born- 
Green equation (eq.[SJ, symbols for the right hand side, eval- 
uated using the full triplet distribution function g^ (crosses) 
and the distribution function g^ (open circles) which is based 
on Kirkwood's superposition approximation. 



volume we consider the Born Green equation [Fj 



3w (2) (ri 2 ) du{r 12 ) 



dri 



dri 



= P 



du(r 13 ) g^(r 1 ,r 2 ,r 3 ) 
dri g{r 12 ) 



dr 3 



(2) 

relating the difference between the mean force and the di- 
rect pair-force to an integral over the force on particle 1 
due to a third particle at r 3 , weighted by the probability 
pg^ dr 3 j 'g(r 12) of finding this particle in dr 3 at r 3 when 
it is known that other particles are located at r\ and r 2 ■ 
This equation is exact if pairwise interactions can be as- 
sumed. To illustrate the importance of three-particle cor- 
relations, we numerically computed the right-hand side 
of eq. J2J using both the full and the approximated triplet 
function, g^ and g*g\, of the T = 4 and T = 46 measure- 
ment and compared it in Fig. (JjJ to the left-hand side of 
eq. (J2J, evaluated using u(r) and g(r). For the strongly 
interacting system (r = 46) the KSA fails completely. 
Three-particle correlations are thus seen to be important 
not only to obtain locally the correct structure, but also 
to obtain globally the correct difference between mean 
and direct force via the Born-Green equation. 
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